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Abstract 



We estimate the minimal mass of selfgravitating polytropic disks using 

the famous Sobolev inequality. This bound resembles the well known mass 

formula for Lane-Emden stars. For ideal gas with the polytropic index n — 3 

the minimal mass is not smaller than the Jeans mass. The accuracy of the 

Oh estimate is verified in a number of numerical examples. The bound works 

C/j well for heavy selfgravitating disks and is less useful for light disks. 
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1 Introduction 

b 

Models of static spherically symmetric configurations of selfgravitating polytropic 
fluids can be reduced to the analysis of the so-called Lane-Emden equations. They 
fall into the class of nonlinear ordinary differential equations. Their solutions can 
be obtained by means of simple numerical methods or even analytically for some 



£S) polytropic exponents. 

lO In contrast to that, equations describing axially symmetric models of selfgrav- 

itating polytropic fluids undergoing a stationary rotation are in general intractable 
analytically. They can be written in the form of the Poisson equation for the grav- 
itational potential and equations of hydrodynamics that can be formally integrated 
to yield an algebraic relation between the specific enthalpy, angular velocity, and 
the gravitational potential. Another possibility, less useful computationally, but 
revealing the mathematical structure of the problem, is to express them as a sin- 

S^ gle nonlinear elliptic equation for the specific enthalpy (or the mass density) with 

Dirichlet conditions imposed on the a priori unknown boundary of the disk. 

We derive analytical results on accretion disks using simple physical informa- 
tion and certain functional analytic methods. There exists the so-called Sobolev 
inequality [1]. In this paper we demonstrate that it can be used in order to estimate 
the mass of selfgravitating toroids. This approach works for those rotation laws 
that possess a centrifugal potential O c such that AO c < (here A is the Laplacian). 
The equation of state of fluid is p - Kp l+l ^ n , where p is the pressure, p the density, 
and n > 3. In analogy to the well known Lane-Emden case, the estimate requires 



the knowledge of the maximal density and temperature of the configuration. The 
same bound is obtained for rotating polytropic stars. It yields rigorously, for ideal 
gas with the polytropic index n = 3, that the mass of stationary systems is not 
smaller than the Jeans mass. That is probably the first rigorous derivation of the 
Jeans inequality. 

The order of the paper is as follows. In the next section all relevant equations 
are displayed. Section 3 gives the estimate for the mass of a selfgravitating disk. In 
section 4 we discuss the spherically symmetric case. There emerges a striking sim- 
ilarity between the expression for the mass of a Lane-Emden star and the estimate 
that is proved in Section 3. Section 5 compares results of the numerical solutions 
describing axially symmetric disks and the analytical estimate of this paper. This 
comparison shows that our approach is particularly robust in the strongly nonlinear 
regime, when the disk masses are much larger than the central mass. Final Section 
contains a brief summary. 

2 Notation and equations 

Consider a disk of perfect fluid rotating around a central point mass. In this case 
stationary Euler equations can be written as 

(U • V)U = -VO - -Vp. (1) 

P 

Here U denotes the fluid velocity; O is the gravitational potential; p denotes the 
density and p the pressure of the fluid. For a selfgravitating disk the gravitational 
potential satisfies 

AO = 4nGp, (2) 

where G is a gravitational constant. 

Let (r, (p, z) denote cylindrical coordinates. We will consider purely rotating, 
axially symmetric configurations, so that U = cod^, where co is the angular velocity. 

In order to find solutions of the above system of equations it is customary to 
assume a fixed form of the rotation law a) = to(r) and the equation of state p = p(r). 
There is a vast literature on the numerical solutions of such problem (see e.g. [2]). 
In this paper we are interested in establishing general analytic bounds on the mass 
of the disk. 

In what follows we will specialize to the polytropic equation of state p = 
Kp l+l / n , where K and n are constant. Let us introduce the specific enthalpy h, 
so that dh = dp I p. For the polytropic equation of state h = K{\ + n)p l ^ n . 

Computing the divergence of Eq. ([TJ yields 

Ah = -AnGp + -d r (rV) - -Ch" - AO c , (3) 

where we have introduced the centrifugal potential 

O c = - f dr'r'o?{r'), (4) 



and a constant C = 4nG/(K(l + ri)) n . 

Notice that Eq. Q can be further simplified by assuming the so-called v-const 
rotation law, i.e, a> = vo/r, where vo is a constant. In this case the centrifugal term 
on the right-hand side vanishes, and we have 

Ah = -Ch n . (5) 

This simple form does not imply that the rotation does not influence the structure of 
the disk. It only means that <1> C is a harmonic function inside the volume occupied 
by a disk. Let us point out that Eq. ([5]) is still difficult to solve, because the boundary 
condition h = is to be posed on an unknown boundary of the disk. 

In this paper we will reserve the symbol Q. for the domain in R. occupied by 
the disk. The disk boundary will be denoted by dQ.. 

Eq. (|5]> is also valid for a static polytropic star. In that case, if we assume 
isotropy, introduce spherical coordinates and properly rescale variables, it reduces 
to the well known Lane-Emden equation. 

3 Estimates of the disk mass 

Assume that A<l> c < 0. For the class of rotation laws of the form a) = const/r p this 
implies that p < 1, so that the v-const rotation is a limiting case. 

Let us multiply both sides of Eq. Q by h and integrate over Q.. It is easy to 
observe that 

- \ d 3 xhAh = | d 3 x\Vh\ 2 = C \ d 3 xh" +l + \ d 3 xhA® c 

JQ. Jfi Jfi JO. 

< C f d 3 xh n+l , (6) 

Jn 

where the left-hand side has been integrated by parts, and we have used the fact 
that h = on d£l. Further steps are adapted from [3]. The last integral in Eq. ([6]) 
can be estimated making use of the Holder inequality. For n > 1 we have 



f d 3 xh n+1 = f d 3 xh"- 1 h 2 < IIV-MLa,™ \\h 

Finally with the help of the Sobolev inequality 

\\h\\ L6(a) <C(3,2)Wh\\ L 2 (a) 



Il 3 (H) 



we arrive at 



Thus 



liv/zii^^c^o.^ll^^l^iiv^ii^. 



\h^l^,Md 3 xh^ > l 



lL 3 / 2 (n) 



Q / -CC 2 (3,2)' 



The Sobolev constant C(3, 2) = 4^ 3 / ( ^/3n 2 ^ 3 ) is a universal number in K. 3 (cf. H] 
4J). The specific enthalpy h can be extended to a function defined on K 3 by setting 
h - outside Q. Such an extension belongs to W Q ' (K 3 ), i.e., the closure of the set 
of compactly-supported C°°(R 3 ) functions in the Sobolev space W 1,2 (R 3 ). 
The mass of a disk is given by 



Let / be 



M = \ d 3 xp = — \ d 3 xh n . 
J n AnG J n 

/ = d 3 xh— = d 3 xh—h n . 
Ja Jo. 



For n > 3 the value of / can be estimated as 

I < 



h" 2 I d 3 xh n 

"max I " ■*■"■ ! 

Jn 



where /z max denotes the maximum value of the enthalpy within the disk. The re- 
versed inequality holds for n < 3. 

A combination of the above results gives the lower bound on the mass of the 
disk in the form 

/ i— «z3 \ _1 I K(\ + n)\ 2 "-3 

M>(47rGVcC 3 (3,2)Vax) = I \ nG 1 P m ^'C- 3 (3,2) (7) 

valid for n > 3 and A<i> c < 0. Here, similarly, p max denotes the maximum of the 
density within the disk. 

The obtained result can be also understood as a bound for the maximum tem- 
perature in the gas configuration. For the ideal gas T = p/um p / ipk^), where // 
is the mean molecular weight, m p denotes the mass of the proton, and &b is the 
Boltzmann constant. Inequality (|7]) can be now written as 

T max < 47rG/im p M3pi ax C 2 (3,2)/ ((1 + n)k B ) , 

where T max denotes the maximal temperature of gas. 

Yet another consequence of ([7]) can be obtained for the value n = 3 of the 
polytropic index. Let p and f denote volume averaged mass density and tempera- 
ture, respectively. For the ideal gas one gets, after simple calculations involving a 
Holder inequality: 

T<K^pi.. (8) 

Inserting that into ([7]) yields the estimate 

3 _ 3 

M>(-^) 2 ^C- 3 (3,2), (9) 
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Figure 1: The integral M as a function of n. The horizontal line represents the 
value of C" 3 (3, 2) = 12.821. 



or, writing C (3,2) explicitely, 
M > 



k B 



4 \Gfim p J yfp 



Ifl 



= M S . 



(10) 



The Jeans mass is usually expressed as 



Mj=j3 



G/UMr 



2 Jj 



(11) 



where /3 is a constant (dependent on the convention of the Jeans mass) of the order 
of 1/2. Bound systems should possess a mass M > Mj, according to a derivation 
that traditionally suffers from severe gaps. Notice that the Jeans mass Mj is close 
to M s , the right hand side of ( 10 1; indeed, M s = Mj x 3 V3~7r/(4/3). Therefore for 
bound systems we obtain rigorously M > 3 V37r/(4/3)Mj. 



4 Lane-Emden stars 

Inequality ([7]) is obviously satisfied for static polytropic stars, described entirely in 
terms of the so-called Lane-Emden functions. 

Assuming the spherical symmetry Eq. Q can be reduced to the Lane-Emden 
equation 



hmv^- 



(12) 



where 

* \ JT(l+n) r ' 

P - PmaxO", and r denotes the distance from the center of the star. Since p max = 
p(r - 0), we have 9(t; = 0) = 1. A textbook exposition of the theory of Lane- 
Emden equation can be found in [5]. 

Let the first zero of 9 (if present) be denoted by £o- If the Lane-Emden function 
corresponding to a given index n has no zeros, we assume £o - °°- A radius 
corresponding to £o will be denoted by R. The mass of the star can be computed as 



3 

(K(\ + n)\ 2 »-3 
\ nG PnZl M, (13) 



Jo 



where 

M = J%^V^ = -4^ 2 |) . 

Since analytical expressions for 9 are only known for n = 0, 1 and 5, the values of 
M have to be computed numerically. 

The similarity between Eqs. and ( 13 1 is remarkable. We see that for n > 3 



there must be M > C 3 (3, 2). 

The general bound on M can be obtained by observing that M attains its min- 
imum value for n = 0. In this case the solution 9 is known analytically and 
M = M m m - 16 V67r/5 = 24.625. Numerical values of M for different indices 
n are shown on Fig.[lJ The factor M is close to C~ 3 (3, 2) for n close to zero, which 
is outside the validity zone of our estimate. The solutions of the Lane-Emden 
equations are, however, relatively easy to be obtained numerically, and this case 
can only be treated as a demonstration of the validity of Eq. (|7]). In the next section 
we will present some numerical results for the selfgravitating disks. 

5 Selfgravitating disks 

In the present Section we shall test the quality of the estimate (|7]) by comparing it 
with appropriate numerical solutions describing selfgravitating rotating disks. The 
numerical method used here is quite standard (2j|6j|7]|. Here we will only sketch 
the main idea. 

Eqs. ([TJ can be formally integrated yielding 

h + <D + <D C = C, (14) 

where <1> C is given by Q and cD satisfies (|2]). The integration constant C is im- 
portant, and its value depends on the solution. Gravitational potential O can be 
expressed in terms of the Green function for the Laplace operator 



X 



GM C C , , p(x') 

<D(x) = —^ - G d 3 x'P-^- (15) 

x - x'\ 
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Figure 2: An example of the density distribution in the selfgravitating rotating disk. 
The plot shows a cross section of the upper hemisphere in a meridian plane. Spatial 
dimensions are expressed in solar radii /?©. The density is color coded in g • cm -3 . 



where M c denotes the central mass. In the implementation of the numerical method 
the singularity of the Green function at x' = x causes a problem, that we avoid in 
the standard way — by expanding the above integral in terms of Legendre polyno- 
mials 0. 

We assume that the domain Q. ranges from r = r m to r - r out in the equatorial 
plane, and the maximal density within the disk reaches a fixed value p max - The 
rotation law is co(r) = vq/t, and the equation of state has the form p = Kp [+[ ^ n . 

The choice of the initial parameters is to a large extent arbitrary. We fix values 
of r m , r ou t, n, p max and M c , whereas constants C, vo and K are computed once the 
solution is known. For vq we have 



V n = 



ln(/outMn) 



where we use the fact that the enthalpy h vanishes on dQ.. Values of C and K can 
be obtained form Eq. ( [14] ) taken at (r m ,z = 0) (or (r out ,z = 0)) and the point where 
P = Pmax respectively. 

The structure of the disk is obtained by an iterative procedure. One starts by 
assuming an initial density distribution in the toroidal shape ranging from r m to 
r max with the a priori given maximum value p max - Then the gravitational potential 



Table 1: Disk masses and their estimates computed from (|7]). Configurations listed 
in the table were computed for n = 3 and M c = IMq. 
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[Rq] Pmax [g ■cm j ] M [Mq] Mass estimate [Af Q ] 
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10 


0.1 


1.4 


5.1 • 10~ 2 


1 


10 
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25 


2.7 


50 


100 


lO" 2 


1.8- 10 2 


7.3 


500 


1000 


10" 5 


1.8- 10 2 


7.3 


2500 


5000 


10~ 6 


2.2 • 10 3 


92 


1 


1000 


10" 4 


1.5- 10 3 


3.8 • 10 2 


1 


100 


0.1 


1.8- 10 3 


3.8 • 10 2 


1 


5000 


10" 5 


1.8- 10 4 


5.5 • 10 3 



is found from the formula ( 15 ). In the next step we compute constants vq, C and 



K corresponding to the assumed density distribution. Finally Eq. ( 14 1 can be used 
to obtain a new approximate solution for h and the corresponding distribution of p. 
If Eq. ( p"4| ) gives a negative value of h in some region of interest then we assume 
that h = there. In this way a new shape of the disk is obtained. This three-step 
procedure is iterated until a satisfactory convergence is reached. 

The disadvantage of such approach is that a good spatial resolution is required 
in order to obtain accurate solutions, and computations of the gravitational poten- 



tial O given by Eq. ( 15 ) become time consuming. 

An example of a disk solution obtained by the procedure described above is 
shown on Fig. [2] The solution was obtained for r m = IRq, r out = 100/?q> n - 3, 
Pmax - O.lg • cm" 3 and M c = IMq. The corresponding mass of the disk is equal 
to M = 1.8 • 10 3 %, and the bound given by Q is M > 3.8 • 10 2 Rq. Here Mq and 
Rq denote the solar mass and radius respectively. 

Table [T] summarizes results obtained for a couple of numerical solutions and 
analytic estimates. The first three columns show values of the initial parameters: 
the inner radius r m , the outer radius r out and the maximal mass density p max , re- 
spectively. The last two columns display the mass of a disk and its estimate derived 
from ([T]). Let us point out that the geometric data, r m and r out , are needed only for 
the numerical calculation. 

The obtained inequality is never saturated. Relatively heavy disks, with the 
mass exceeding the central mass by 3-4 orders of magnitude, have masses close to 
that predicted by (|7]). It is clear that the accuracy of the functional-analytic bound 
increases with the mass of the disk, and thus with the increase of the selfgravity. 
The exact numerical mass and the predicted mass differ by a factor of three in the 
heavy end of masses, and by a less than two orders of magnitude for light disks. 



6 Conclusions 

We derive, using certain functional inequalities, an analytic estimate for the mass 
of selfgravitating axially symmetric stationary configurations of polytropic fluids. 
It is valid for polytropic indices n > 3, both for rotating stars and for accretion 
disks with centrifugal potentials satisfying the condition AO c < 0. This class of 
potentials is quite general and includes several common types of rotation, with the 
rigid co = const and the v-const rotations. 

The bound on the mass is given in terms of the polytropic index, maximal 
density and maximal temperature of the gas. The obtained expression is strikingly 
similar to the mass formula of the Lane-Emden stars. One can use this result 
in order to obtain the Jeans inequality for equilibrium ideal gas disks and non- 
spherical stars for the polytropic index n = 3. We believe that this is the first 
rigorous derivation of the Jeans mass and of the Jeans inequality for stationary 
systems. 

The accuracy of the analytic bound increases with the mass of the disk. The 
exact numerical mass and the predicted mass differ by a factor of three in the heavy 
end of disk masses, and by a less than two orders of magnitude for lighter disks. 
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